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Abstract 

Large-scale research programmes seeking to characterize the C 4 pathway have a requirement for a simple, high 
throughput screen that quantifies photorespiratory activity in C 3 and C 4 model systems. At present, approaches rely 
on model-fitting to assimilatory responses {A/C, curves, PSII quantum yield) or real-time carbon isotope discrimina- 
tion, which are complicated and time-consuming. Here we present a method, and the associated theory, to determine 
the effectiveness of the C 4 carboxylation, carbon concentration mechanism (CCM) by assessing the responsiveness 
of V 0 IV C , the ratio of RuBisCO oxygenase to carboxylase activity, upon transfer to low 0 2 . This determination com- 
pares concurrent gas exchange and pulse-modulated chlorophyll fluorescence under ambient and low 0 2 , using 
widely available equipment. Run time for the procedure can take as little as 6 minutes if plants are pre-adapted. The 
responsiveness of V Q IV C is derived for typical C 3 (tobacco, rice, wheat) and C 4 (maize, Miscanthus, cleome) plants, and 
compared with full C 3 and C 4 model systems. We also undertake sensitivity analyses to determine the impact of R L ight 
(respiration in the light) and the effectiveness of the light saturating pulse used by fluorescence systems. The results 
show that the method can readily resolve variations in photorespiratory activity between C 3 and C 4 plants and could 
be used to rapidly screen large numbers of mutants or transformants in high throughput studies. 

Keywords: C 4 , C 3 , photosynthesis, RuBisCO, oxygenation, carboxylation, carbon concentration mechanism (CCM), Cleome 
gynandra, rice, maize, wheat, Miscanthus. 



Introduction 

In most photosynthetic organisms Ribulose Bisphosphate 
Carboxylase Oxygenase (RuBisCO) catalyses the first key 
step in carbon assimilation, reacting ribulose- 1,5-bisphos- 
phate with C0 2 to produce two molecules of 3-phosphoglyc- 
erate (PGA). Oxygen competitively inhibits this reaction and 
leads to the synthesis of the 2-carbon compound phosphogly- 
collate, which is recycled to PGA (consuming ATP, and then 
NADPH) and C0 2 by the photorespiratory cycle (Yoshimura 
et al., 2004; Sage et al, 2012). The result of photorespiration 
is a noticeable carbon loss and a consequent metabolic cost 
for carbon recapture and for the recycling of photorespira- 
tory intermediates (Ehleringer and Pearcy, 1983; Pearcy and 
Ehleringer, 1984; Eckardt, 2005). Many plants have evolved 
strategies to reduce photorespiration by increasing the level 
of C0 2 around RuBisCO, including both crassulacean acid 



metabolism (CAM) and the C 4 photosynthetic pathway 
(Dodd et al. , 2002; Sage, 2004; Sage et al. , 201 1 ; Osborne and 
Sack, 2012; Griffiths et al, 2013; Owen and Griffiths, 2013). 
C 4 photosynthesis is most often based on a two-celled carbon 
concentrating mechanism, where HC0 3 is first fixed into the 
four-carbon compound oxaloacetic acid (OAA) in the meso- 
phyll by phosphoeno/pyruvate carboxylase (PEPC). OAA is 
then reduced to malate or transaminated to aspartate and the 
resulting C 4 -(amino)acid is shuttled into the bundle sheath 
(BS), where it is decarboxylated, releasing C0 2 for refla- 
tion by RuBisCO (Hibberd and Covshoff, 2010; Bellasio and 
Griffiths, 2014c). 

Although the enzymes catalysing the core C 4 carbon con- 
centration mechanism (CCM) are well characterized (Kanai 
and Edwards, 1999), many of the genes responsible for the 
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accompanying anatomical alterations or for generating and 
maintaining expression of the C 4 cycle genes (Hibberd et al, 
2008; Langdale, 201 1) have yet to be identified. One approach 
that is increasingly proving useful to identify candidate genes 
underlying the C 4 pathway is comparative transcriptom- 
ics of samples either undergoing C 3 or C 4 photosynthe- 
sis (Brautigam et al, 2011; Gowik et al, 2011; John et al, 
2014), or tissues in the process of inducing the full C 4 system 
(Li et al, 2010; Pick et al, 2011; Chang et al, 2012; Wang 
et al, 2013). Because stable transformation of C 4 species is 
typically time-consuming, introduction of RNA interference 
constructs via a transient Agrobacterium tumefaciens-ba.sed 
system would be very helpful in screening these candidates 
being generated from transcriptomics. 

At present, techniques used to screen for mutants possess- 
ing defective, or enhanced CCM characteristics are time-con- 
suming (Table 1). Analysing the response of assimilation (A) to 
decreasing C0 2 concentration in the substomatal cavity (Q), 
as A/Q curves (Long and Bernacchi, 2003; Yin et al, 2011a) 
can take 45 minutes per replicate leaf, and an appropriate 
model, which may require a priori knowledge of species-spe- 
cific limitations (Laisk and Edwards, 2000; von Caemmerer, 
2000, 2013; Yin and Struik, 2009; Yin et al, 2009; Yin et al, 
2011b). 13 C/ 12 C discrimination during photosynthesis (Evans 
et al , 1986) can also be used, and a comparison with stomatal 



conductance allows the internal mesophyll conductance, or 
extent of CCM or PEPC activity, to be resolved (Meyer et al, 
2008; Kromdijk et al. , 2010; Pengelly et al. ,2010; Bellasio and 
Griffiths, 2014a, b, c). However, this latter technique is sensi- 
tive, and requires either off-line sample preparation for mass 
spectrometric analyses or specialized laser equipment which 
is not readily available (Table 1). 

In this paper we describe a novel method, and present the 
associated theory, to determine rates of photorespiration 
from instantaneous rates of RuBisCO carboxylation and oxy- 
genation. The approach compares concurrent gas exchange 
and pulse-modulated chlorophyll fluorescence measurements 
under ambient and low 0 2 . Under these non-photorespira- 
tory conditions assimilation (A) increases, because RuBisCO 
competitive inhibition from 0 2 is reduced. In contrast, Y( II) 
decreases because the demand for NADPH associated with 
photorespiratory by-product cycling (and reduction) is lower, 
and cannot entirely be offset by the increase in A. The new 
method combines developments in approaches using gas 
exchange (Sharkey, 1988; Long and Bernacchi, 2003; Ripley 
et al , 2007) and the quantitative interpretation of quantum 
yield (Yin et al, 2004, 2009, 2011b; Yin and Struik, 2009, 
2012; Bellasio and Griffiths, 2014b). This new method can be 
performed with off-the-shelf commercial equipment, which 
is generally available in ecophysiology laboratories. The 



Table 1 . Comparison between methods screening for activity of a functional CCM 



Method 


Advantages and limitations 


Reference 


Dry matter isotopic discrimination 


'Specialized equipment 

Integrates the isotopic signal throughout 

growth 

*Cannot resolve transient changes in 
assimilatory physiology 


Cernusak era/. (2013) 


On line isotopic discrimination 


'Laser is no longer commercially available 


Evans ef al. (1 986); Bellasio and 




'Maintenance costs of isotope ratio 


Griffiths (2014b); von Caemmerer 




mass spectrometer 


ef al. (201 4) 




*Need of highly skilled operator 






'Difficult computation and parameterization 




A/Q curves 


•Requires a priori knowledge of the 


Long and Bernacchi (2003); 




limitations underpinning each part for the 


Yin ef al. (2009) 




A/Cj curve for correct model fitting 






'Result may depend on experimental 






routine 




Gas exchange and fluorescence 


'Requires initial response curve for 


Long and Bernacchi (2003); 




parameterisation 


Martins ef al. (2013) 




'Requires model fitting 




0 2 sensitivity of carboxylation efficiency 


'Delicate experimental routine 


Laisk ef al. (2002); Yin ef al. (2009) 


Assimilation increase under low 0 2 


'Ease of determination 

'Ignores the effect of changing 0 2 

concentration on Y(ll) 


Sharkey (1988); Ripley ef al. (2007) 


Gas exchange and fluorescence 


'Rapid (6 minutes) 

'Widely available equipment 

'Independent of leaf size 

'Ease of determination and calculation 

'Does not require fitting or 

parameterisation 

'Assessment under growth conditions 


This study 
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procedure takes as little as 6 minutes to perform if plants are 
pre-adapted, making it significantly faster than A/Q curves 
and potentially useful as a high-throughput approach for 
assessing C 4 activity in mutant screens, the progeny from C 3 
C 4 crosses or C 3 -C 4 intermediates. 

Materials and methods 

Plants 

Plants of Miscanthus (Miscanthus giganteus), cleome (Cleome gynan- 
dra), maize (Zea mays L.), wheat (Triticum aestivum L.), tobacco 
{Nicotiana tabacum L.), and rice (Oryza sativa L.) were grown at 
the Plant Growth Facility located at the University of Cambridge 
Botanic Garden in controlled environment growth rooms (Conviron 
Ltd, Winnipeg, Canada) set at 16h day length, temperature of 
25 °C/23 °C (day/night), 40% relative humidity, and photosynthetic 
photon flux density (PPFD)=300 umol mr 2 s -1 . Plants were manu- 
ally watered daily, with particular care to avoid overwatering. 

Gas exchange measurements with concurrent PSII yield 

Measurements were performed with an infra-red gas analyser 
(IRGA, a LI6400XT, Ll-cor, USA), fitted with a 6400-40 leaf 
chamber fluorometer. The IRGA was fed with C0 2 (through 
the IRGA gas mixing unit) and ambient air. Gas flow was set at 
1 50 umol s~' . Reference C0 2 was set at 200 umol moL 1 (Figure 1 
and Table 1) or set alternatively at 400, 300, 200, 150, 100, and 50 
umol s -1 (Figure 3). Block temperature was controlled at 35 °C. 
The fluorometer was set to multiphase pulse with factory setting, 
target intensity=10 and ramp depth=40% (Loriaux et at, 2013). 
A portion of a light-adapted leaf was clamped in the cuvette. The 
leaf was allowed to reach stable photosynthetic conditions under 
PPFD=300 umol mr 2 s _1 (factory setting: 90% red, 10% blue). 
Photosynthesis was measured every 10 s for 30 s (the three values 
were then averaged) and a multiphase pulse was applied for the 
determination of Y(II). A humidified 2% 0 2 /N 2 gas (pre-mixed, 
BOC, Guilford, UK) was switched to supply the inlet of the IRGA. 
The gas was allowed to completely flush the cuvette (c. 6min). 
Photosynthesis was measured every 10 s for 30 s (the three values 
were then averaged) and a multiphase pulse was applied for the 
determination of Y(II). Light was turned off, the inlet was fed 
with ambient air, the reference C0 2 was set at 500 umol mol , 
similar to the lab C0 2 concentration (c. 550 umol mol -1 ) to mini- 
mize the errors caused by C0 2 leakage (Boesgaard et cil., 2013), 
and flow was set to 40 umol s -1 . Once the cuvette had been flushed, 
and the signal stabilised (c. 5min), respiration was measured every 
10 s for 2min (the values were then averaged). C a was not adjusted 
to account for changes in stomatal conductance or for the control 
of Q during this procedure. This avoided the need for IRGA rec- 
alibration as the Y(II) measurements are independent of Q. The 
measured A and Y( II) under low and ambient 0 2 , together with an 
estimate of Flight (see below), were used to determine RuBisCO 
rate of carboxylation (V c ), RuBisCO rate of oxygenation (Vq), 
and the rate of photorespiratory C0 2 evolution in the light (F). 

Theory 

RuBisCO catalyses two reactions: a carboxylase reaction whereby 
Ribulose BisPhosphate (RuBP) is carboxylated to form two mol- 
ecules of phosphoglyceric acid (PGA), and an oxygenase reaction 
whereby RuBP is oxygenated to form one PGA and one glycol- 
late molecule. Each carboxylase event requires 2 NADPH for the 
reduction of the 2 PGA molecules formed. Each oxygenase event 
requires 1 NADPH for the reduction of the PGA directly produced 
by RuBisCO, 0.5 NADPH to recycle glycollate, and 0.5 NADPH 
to reduce the PGA regenerated, which total 2 NADPH (Bellasio 
and Griffiths, 2014c). The overall NADPH demand, at steady-state, 



equals the total photosynthetic NADPH production rate / N adph 
(Yin et aL, 2004; Yin and Struik, 2012): 



= 2K C +2K 0 



(1) 



Where /nadph is the total NADPH produced for photosynthesis, 
V c is RuBisCO carboxylation rate, and V Q is RuBisCO oxygenation 
rate. Notably, this reducing power requirement is the same for all 
types of photosynthesis, as active types of CCM require additional 
ATP but not NADPH. In line with von Caemmerer (2000) equa- 
tion 1 assumes that PGA is entirely reduced, and therefore the small 

quantity of PGA consumed by respiration (- Flight) is neglected, 

in fact under growth light irradiance 2 V c +2 V Q »^ Flight* unless at 

very low irradiances, see equation 7 in Bellasio and Griffiths (2014c). 

Although the carboxylation reaction of RuBisCO consumes C0 2 , 
the regeneration of glycollate releases 0.5 C0 2 for each oxygenase 
catalytic event. C0 2 is also produced by light respiration, a process 
which is active during photosynthesis to support basal metabo- 
lism. The net assimilation rate (A, which is the quantity measured 
through gas exchange) results from summing the C0 2 consumed by 
RuBisCO, the C0 2 produced by glycollate regeneration and the C0 2 
produced by respiration: 



A = V r 



(2) 



Where A is net C0 2 assimilation, Flight is respiration in the light 
and other variables were previously defined. Notably, this equation 
is universal for all types of photosynthesis (von Caemmerer, 2013). 

For the definition of gross assimilation {GA-A + ^ L ight)> equa- 
tion 2 can be rearranged: 



V c =GA + \v Q 

Equation 1 and 3 can be combined to give: 

T7 _ /nadph ~ 2GA 
3 



(3) 



(4) 



The rate of photorespiratory C0 2 evolution, F can be calculated 
as: (von Caemmerer, 2013) 



(5) 



Under low 0 2 , V Q can be approximated to =0, hence, from 
equation 4: 



J 



NADPH Low 02 



■2GA, 



(6) 



Which is valid when V o ~0. 

NADPH is produced through linear electron flow. Independently 
from where this reaction is located (e.g. in mesophyll cells), elec- 
trons are invariably extracted from water by PSII (Yin and Struik, 
2012), therefore /nadph is proportional to Y(II) (Yin and Struik, 
2012). This allows /nadph to be calculated under photorespiratory 
conditions using the information derived under non-photorespi- 
ratory conditions, and can be expressed as (Bellasio and Griffiths, 
2014b): 



J* 



, = /, 



NADPH Low 02 



Y(JI) 



(7) 



Where /nadph and Y( II) refer to ambient 0 2 conditions. Equation 
7 has been validated in C 3 and C 4 plants (Yin et al, 2009, 2011b; 
Bellasio and Griffiths, 2014b, c) but it is worth noting that equation 
7 is a mathematical simplification and holds true when: (i) photores- 
piration is negligible under non-photorespiratory conditions, which 
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is a widely used simplification; (ii) -R LIGH t does not vary between 
low and ambient 0 2 — this is also a fair assumption because any 0 2 
effect is generally negligible (Badger, 1985; Gupta et al., 2009); (iii) 
the allocation to alternative sinks (non-assimilatory and non-pho- 
torespiratory) is proportional to Y(II). This is the normal case in C 4 
plants where the relationship between Y(II) and Y( C0 2 ) has a null 
intercept (Edwards and Baker, 1993). When that is not the case, for 
instance when the allocation to alternative sinks is constant, equa- 
tion 7 would also hold true if the allocation to alternative sinks is 
small compared with Y(II). This is the normal case in C 3 plants 
(Valentini et al, 1995; Martins et al, 2013). Should the allocation 
to alternative sinks be large, equation 7 would still hold true math- 



ematically when 



Y(U) 



is close to the unity. The implications 



o: 

for method accuracy are detailed in the discussion. 
Equation 3, 4, 6, and 7 can be combined to obtain: 



2GA Y(II) 
Vo_ LOW ° 2 Y(II) Lm 



-2GA 



(8) 



Y(H\ 



Which expresses the RuBisCO rate of oxygenation relative to 
carboxylation. The influence on the quality of Flight estimate on 
V 0 /V c is described in the discussion, together with the other factors 
influencing the results. 



Modelling C 3 and C 4 Vc/V c 

The data obtained for tobacco and maize were compared with a 
simulated V Q IV C based on the validated von Caemmerer models 
for C 3 and C 4 photosynthesis. Briefly, for tobacco, the response of 
A to Q was modelled using the quadratic equation (Table 3, equa- 
tion 9) proposed by Ethier and Livingston (2004), which takes into 
account mesophyll conductance to C0 2 . The C0 2 concentration at 
the site of carboxylation C c was then calculated through the supply 
function of mesophyll (equation 10), and, finally V Q /V C was simu- 
lated from the kinetic properties of RuBisCO and the ratio between 
C c and the 0 2 concentration at the site of carboxylation (equation 
11). For maize (Table 4), firstly we simulated the responses of V P 
and A to decreasing Q, using the equations for the enzyme-limited 
model for C 4 photosynthesis (equation 12 and 16, respectively). 
These were used to simulate the C0 2 and 0 2 concentration in the 
bundle sheath (equation 13 and 14, respectively), the ratio of which, 
together with RuBisCO specificity, was used to simulate V Q IV C 
(equation 15 and 17). 

Results 

Figure 1 displays a typical primary data profile for a C 3 
tobacco leaf, showing the interaction between steady state 
assimilation (A) and quantum yield of PSII, Y(II), during 
the transition from ambient to low 0 2 (21 to 2% 0 2 ), with 
hatched areas indicating the steady state conditions under 
which readings were taken to derive V Q /V C . Under non- 
photorespiratory conditions, A increases because of the 
lower competitive inhibition of 0 2 , whereas Y( II) decreases 
owing to the lower NADPH demand for photorespiratory 
by-product recycling and reduction. The experimental con- 
ditions were deliberately chosen to minimize reductions of 
quantum yield at saturating light (relatively low PPFD of 
300 umol nT 2 s -1 ), and enhance photorespiratory responses 
to low 0 2 partial pressure (measurements at 200 umol 
moF 1 C0 2 ) (Fig. 1 and Table 2). Subsequently, V Q /V C was 




M 3.0 



100 200 300 400 500 600 
Time / s 

Fig. 1. Summary of experimental approach. One representative dataset 
from C 3 tobacco is presented. Once stable assimilatory conditions 
are reached, a first set of data are recorded (left hatched area). The 
background gas is then switched from ambient to 2% 0 2 . After a suitable 
acclimation time to allow flushing of the cuvette and reacclimation (c. 
6min), a second set of data are recorded (right hatched area). The 
response of assimilation (triangles) and Photosystem II yield Y(ll) (squares) 
during the experiment are shown. 



measured on C 3 tobacco and C 4 maize using different C0 2 
concentrations in the reference gas: 400, 300, 200, 150, 100, 
and 50 umol moF 1 (Fig. 2) and results were compared with 
simulated values of VqIV c generated with the validated von 
Caemmerer C 3 and C 4 models. To facilitate the comparison, 
data were plotted against the substomatal C0 2 concentra- 
tion Q. As expected, under decreasing Q, VJ V c becomes 
progressively higher in tobacco but it is only marginally 
affected in maize. The measured data track the trend and 
magnitude of the theoretical curves in C 3 , whereas we could 
not capture the theoretical increase in VJ V c expected when 
Q was close to zero. This may be due to errors in the deter- 
mination of Q at very low stomatal conductance or to the 
simplifications used to resolve equation 7. Our data slightly 
underestimate V Q /V C derived using pulsed of 13 C enriched 
C0 2 (Busch et al, 2013), which, however, lay above the 
curve simulated with the von Caemmerer C 3 model (see 
Fig. 2). 

Additional measurements were undertaken with the IRGA, 
including a recalibration procedure to account for the chang- 
ing sensitivity to water vapour pressure after the transition 
to low 0 2 , but stomatal conductance was reduced on average 
by 1% and internal C0 2 concentration, Q, by 3 umol moF 1 
(data not shown). In the subsequent sections, primary data 
for VJ V c determinations using this new method (calculated 
from equation 4) are initially presented for three representa- 
tives of C 3 and C 4 species. We then undertake a systematic 
error analysis of the method, to include the impact of biolog- 
ical and environmental variables. These include physiological 
components (Flight) and Fm', as well as light intensity and 
C0 2 concentration used during experimentation. 

Variability between and within populations 

Table 2 demonstrates that the method clearly discriminates 
between C 4 species, possessing a functional CCM, and C 3 
species with higher rates of photorespiration. VJ V c ranged 
from 0.0435 to 0.0852 for the representative C 4 species, with 
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Table 2. Example of variability within populations and between populations displayed by plants with different pathways of assimilation 

Vq/V c was measured on species (Miscanthus, Cleome gynandra, maize, wheat, tobacco, and rice) under photosynthetic photon flux density 
(PPFD) of 300 jumol rrr 2 s"\ and C a =200 //mol mol" 1 . 



Population 


n 


Mean V 0 /V c 


Standard deviation 


Coefficient of variation 


Miscanthus 


7 


0.0504 


0.0091 


18% 


Cleome gynandra 


5 


0.0852 


0.0046 


5.4% 


Maize 


4 


0.0435 


0.0074 


17% 


Wheat 


3 


0.522 


0.071 


14% 


Tobacco 


4 


0.533 


0.030 


5.5% 


Rice 


4 


0.569 


0.037 


6.5% 





^ Tobacco 




Maize 




C3 Model 




C4 Model 


\ A 




■ \ * 












■ \ "T^- — . 








" 1 — 




50 100 150 200 250 


300 


C j 1 nmol mol 1 





Fig. 2. Vq/V c measured under different C0 2 concentrations in 
the substomatal cavity (Q), obtained by imposing reference C0 2 
concentrations of 400, 300, 200, 1 50, 1 00, and 50 umol mol" 1 for 
C 3 tobacco (triangles) and C 4 maize (squares). Data are compared 
with simulated Vq/V c using the validated von Caemmerer C 3 and C 4 
models (lines, see also Table 3 and 4). With decreasing Q, V 0 /V c gets 
progressively higher in tobacco but it is only marginally affected in maize, 
C0 2 concentration can therefore be used to control the resolution of the 
method. All data shown, n=4. 

coefficients of variation ranging from c. 15% down to 5% in 
C. gynandra (Table 2). For the C 3 species, VJV C ranged from 
0.522 to 0.569, with a low coefficient of variation in tobacco 
and rice around 6% (Table 2). The magnitude of the offset 
between C 3 and C 4 systems, if being used as a rapid screen, 
would allow changes in expression of C 4 characteristics to be 
clearly resolved. Such an approach would then allow more 
detailed characterisation of selected transformants, C 2 , or 
C3-C4 intermediates to be undertaken. 

Accuracy of P. UGHT estimates 

To account for the extent that Flight affected the measure- 
ment of Vq/Vc, a sensitivity analysis was used to determine 
how i? LIG HT influences V Q /V C (Fig. 3). To do so, equation 8 
was calculated for a realistic dataset (-Klight = 1 umol nr 2 
s" 1 , V o IV c -0.2 and Y( II) -0.65) at variable assimilation 
values. Then, test values for V Q /V C were calculated after 
Flight was varied to 2 umol nT 2 s _1 (+100%), 1.5 umol 
nr 2 s" 1 (+50%), 1.2 ixmol nT 2 s" 1 (+20%), 0.8 umol nT 2 s" 1 
(-20%), 0.5 umol nr 2 s" 1 (-50%), 0 umol nr 2 s" 1 (-100%, 
GA-A). The deviation from the set V Q /V C value (0.2) rep- 
resented the effect of errors in the evaluation of Flight on 
VJ V c . Figure 3 shows that VJ V c was relatively insensitive 



to Flight : f° r assimilation rates higher than 4 umol nr 2 
s -1 , Flight values which differed ± 50% resulted in an error 
lower than 4% in relative terms. Flight overestimation 
resulted in a lower error than i? LIG HT underestimation. For 
these reasons there is generally no need for a high quality 
estimate of Flight- 

Accuracy of Fm' measurements 

Equations 7 and 8 require the photochemical yield of PSII, 
Y( II) . This is determined according to the formula of Genty 
(Genty et al, 1989; Maxwell and Johnson, 2000; Kramer 
et al, 2004), whereby Y(II) is calculated as the difference 
between the light-saturated chlorophyll fluorescence signal 
{Fm ') minus the chlorophyll fluorescence signal measured 
during photosynthesis (Fs), expressed as relative to Fm'. Key 
to this technique is achieving full saturation of PSII in the 
determination of Fm' (Earl and Ennahli, 2004; Loriaux et al. , 
2006; Harbinson, 2013; Loriaux et al, 2013). Sub-saturating 
light pulses result in the underestimation of Fm'; however, the 
degree of underestimation depends not only on the saturating 
pulse spectra and intensity, but also on the species, the growth 
light intensity, and the light intensity used during the meas- 
urements (Earl and Ennahli, 2004). 

Here, we show how a given Fm' underestimation influences 
the values for VJ V c (Fig. 4). To do so, equation 8 was set to 
physiologically realistic conditions (-Rlight = 1 umol nr 2 s -1 , 
Fo/F c =0.2, and A=5 umol nr 2 s 1 ), at different Y(II) values. 
Underestimates of Fm' were then introduced by multiplying 
the realistic Fm' value by, successively, 0.99 (-1%), 0.98 (-2%), 
0.97 (-3%), and 0.95 (-5%). The difference between the two 
values represented the effect of Fm' underestimation on V Q /V C . 
Figure 4 shows that V Q /V C was sensitive to Fm' underestima- 
tion; for instance the relative error of VJ V c was c. 20% when 
Y(II) was 0.15 and Fm' was underestimated by 3%. The error 
increased hyperbolically at decreasing Y(II), and increased 
proportionally as the Fm' underestimation was increased. 

Light intensity and C0 2 concentration used for 
experimentation 

High light intensities (e.g. PPFD>1000 umol nT 2 s -1 ) result in 
a low PSII yield, which may potentially amplify the systematic 
error from any fm' underestimation (see above). Similarly, small 
Y(II) could potentially lead to VJV C underestimation when 
the allocation to alternative sinks is significant (see description 
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of equation 7). Further, high light conditions require longer 
timescales to reach stable photosynthetic conditions. On the 
other hand, depending on growth conditions, low light intensi- 
ties (e.g. <100 umol nr 2 s -1 ) might lead to low assimilation rates, 
which could amplify the systematic errors in the estimation of 
Flight (see Fig. 3 and above). For these reasons, intermediate 
light intensities represent the best solution, whereby Y( II) and 
A are both high. For instance, values at the top end of the linear 
region of the light response curve would be ideal. These gener- 
ally correspond to the growth light intensity. 

C0 2 concentration in the cuvette (C a ) can be used to manip- 
ulate photorespiration. Figure 2 shows the measured and 
predicted VqJVq of C 3 and C 4 plants under different C0 2 con- 
centrations. Because of the CCM, VJVc is low in maize, even 
at low Q, whereas in wheat VJ V c increases hyperbolically at 
decreasing Q. This contrasting behaviour allows the resolution 
of the method to be manipulated by changing the C0 2 con- 
centration in the background gas. However, decreasing C0 2 




2 4 6 8 10 12 14 16 18 
A / |imol m"V 1 

Fig. 3. Sensitivity to errors in the determination of R L ight- True values were 
simulated by calculating equation 8 for R L ight=1 urnol m ~ 2 s~ 1 , Vq/V c = 0.2, 
and Y(II)=Q.65 at variable assimilation (A) values. Test values of Vq/V c were 
then calculated by solving equation 8 at different values for Alight 
2 umol rrf 2 s _1 (+100%), 1 .5 umol rrf 2 s" 1 (+50%), 1 .2 |imol rrf 2 s _1 
(+20%), 0.8 umol rrf 2 s _1 (-20%), 0.5 |imol rrf 2 s _1 (-50%), 0 |imol rrf 2 s _1 
(-100%, GA=A). The difference in Vc/V c between the test minus the true 
value was expressed as relative to the true value. 
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Fig. 4. Sensitivity to errors in the determination of Fm'. True values were 
simulated by calculating equation 8 for R L ight=1 u.mol nr 2 s~ 1 , Vc/V C =Q.2 
and A=5 umol rrf 2 s" 1 at different Y(ll) values. Test values of Vo/V c were then 
calculated by solving equation 8 introducing increasing Fm' underestimation: 
-1 , -2, -3, and -5%. The difference in Vq/V c between the test minus the 
true value was expressed as relative to the true value. 



concentration is disadvantageous because: (i) low Q results 
in quenching of PSII yield, which may potentially amplify 
the systematic error determined by Fm'' underestimation (see 
above); at the same time (ii) low Y( II) would amplify the mag- 
nitude of VqIVq underestimation owing to the partitioning of 
Y(II) to alternative sinks (see description of equation 7); (iii) 
under low C a , more time is required to reach stable photosyn- 
thetic conditions, which result in lower throughput; (iv) low 
C a increases the driving force for diffusion from outside of the 
cuvette, which may constitute a potential source of error, espe- 
cially when assimilation is low (Boesgaard et al., 2013). For 
these reasons the optimal C a will depend on the purpose of the 
analysis, and on the desired resolution and speed. 

Discussion 

This method is based upon the difference in net assimila- 
tion (A) and photosystem II yield (Y(IIJ) observed when 
the gas supplied to an actively photosynthesizing leaf 
is switched from ambient 0 2 to low 0 2 . The goal was to 
develop a relatively quick, readily available method, which 
could be used to screen large numbers of transformants, 
C3-C4, C 2 , or photorespiratory refixation variants (Busch 
et al 2013; Oakley et al., 2014) in a given population of 
plants. The data show that the method readily distinguishes 
between V Q IV C for typical C 3 and C 4 plants (Table 2), and, 
given the low coefficients of variation, should detect more 
subtle variations in C 4 repression or activation within a 
screen. It would then be possible to subject plants iden- 
tified in this way to a more detailed, conventional gas 
exchange or stable isotope screen, to identify contributory 
morphological, metabolic or genetic factors. In the sub- 
sequent discussion, we explore the theoretical and practi- 
cal limitations underpinning the accuracy of the method, 
and improvements that could be instituted to enhance the 
outputs, if high sample throughput was not a primary 
limitation. 

Other methods have been proposed to determine the 
contribution of photorespiration in vivo through gas 
exchange measurements. The method proposed by Ripley 
et al. (2007) uses only the increase in assimilation under 
non-photorespiratory conditions, and therefore ignores 
the effect on Y(II), In our work we observed that Y(II) 
is generally influenced by changes in 0 2 concentration 
(Figure 1), even in C 4 plants (see Fig. 2 in Bellasio and 
Griffiths, 2014b); therefore it is important to take into 
account the feedback from assimilation on photosystem 
II yield. Long and Bernacchi (Long and Bernacchi, 2003) 
proposed a comprehensive method to determine the parti- 
tioning of total electron transport rate between photores- 
piratory and assimilatory demand. Their protocol requires 
an initial light or AIQ response so as to fit a linear relation- 
ship between quantum yield for C0 2 fixation Y( C0 2 ) and 
quantum yield of photosystem II, Y( II), 

In comparison, the simple method that we have proposed 
requires no previous parameterization, no curve fitting, and 
no knowledge of the underpinning physiology or biochemi- 
cal constants. It is also independent of leaf area, as when 
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deriving V Q /V C from equation 8, both the numerator and the 
denominator are proportional to leaf area, a huge advantage 
for small or dissected leaves. The likelihood of triose phos- 
phate limitation (Sharkey, 1988) is minimized under the rela- 
tively low light intensities and low Q, which are optimal for 
this protocol. The determination of VJ V c could take as lit- 
tle as c. 6 min, although the complete routine was longer (c. 
40min) as leaves were allowed to acclimate before measure- 
ment of both assimilation and dark respiration. Therefore, 
the run time can be minimized by measuring assimilation 
under growth conditions (e.g. at growth light intensity and 
C0 2 concentration), and either measuring respiration after 
all plants have been collectively dark-adapted, or estimating 
it separately (see below). 

Other factors affecting accuracy of Vc/Vg 
determination 

As shown in Fig. 3, the estimation of Flight is impor- 
tant when calculating gross assimilation using eqn. 8 
(GA-A+R UGm ) at low assimilation rates. Flight can be 
determined with several methods; for instance, by linear 
regression of assimilation (^4) versus irradiance (under very 



low irradiance e.g. < 1 50 umol photons m 2 s _1 ), by linear 
regression of A versus irradiance multiplied by Y( II) [under 
moderate irradiance, e.g. <400 umol photons m" 2 s -1 (Yin 
et al, 2011a)], by non-linear regression [throughout the 
light response curve (Prioul and Chartier, 1977; Dougherty 
et al, 1994)] or assumed to equal dark respiration [e.g. 
(Kromdijk et al, 2010; Ubierna et al, 2013)]. These meth- 
ods do not necessarily yield the same Flight values, and so, 
the degree of similarity between different Flight estimates 
depends on the species and growth conditions. For instance, 
in Cocklebur (Xanthium strumarium L., Asteraceae), Flight 
was significantly different from dark respiration (Tcherkez 
et al, 2008), whereas in maize Flight is generally non- 
significantly different from dark respiration (C. Bellasio, 
unpublished data). The most suitable method to estimate 
Flight should therefore be evaluated on a case-by-case basis 
(see Bellasio and Griffiths, 2014a), and for a uniform popu- 
lation (e.g. one species or set of transformants in a growth 
chamber), Flight could be estimated on a subset of indi- 
viduals, with one of the methods described above. If dark 
respiration is used as a proxy, the quality of the estimate can 
be increased using large chambers and low flow rates. In a 
diverse population, i? LIG HT could be estimated by measuring 



Table 3. Model for C 3 photosynthesis 



Symbol 



Definition/calculation 



Equation 



Values/Units/References 



Net Assimilation A = 
where: 
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b JVcmax Flight ) + C| + Kc(1+ 0_ 
9m Ko' 



C=R L IGHT Q+K c (1 + — ) 

-i/n™„ (c-r) 



Ethier and Livingston (2004) 



C0 2 partial pressure at the site of 
carboxylation C c = Q 



(10) 



ixbar 



9m 

K c 
K 0 
O 

Flight 

^Cmax 



C0 2 concentration in the intercellular 
spaces as calculated by the IRGA. 
Mesophyll conductance to C0 2 
RuBisCO Michaelis-Menten constant 
for C0 2 

RuBisCO Michaelis-Menten constant 
for 0 2 

0 2 partial pressure at the site of 

carboxylation 

Respiration in the light 

Maximum RuBisCO carboxylation rate 

Vq _ O 

K3 ^CmaxKo C c 



(11) 



(xmol mol 1 (Li-cor6400 manual equation 1-18) 

0.25 mol rrf 2 s~ 1 bar 1 (Ethier and Livingston, 2004) 
319.3 |ibar (Ethier and Livingston, 2004) 

277 1 00 |ibar (Ethier and Livingston, 2004) 

200000 |ibar 

0.63 umol nr 2 s~ 1 

34.7 umol nr 2 s~ 1 (Ethier and Livingston, 2004) 
equation 2.16 in (von Caemmerer, 2000) 



^Omax 

r* 



Maximum RuBisCO oxygenation rate 
C0 2 compensation point in absence of 
dark respiration 



13.25 umol nr 2 s~ 1 (Ethier and Livingston, 2004) 
44 Libar 
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Table 4. Model for C 4 photosynthesis 



Symbol Definition/calculation 



Equation Values/Units/References 



Net Assimilation A- 
where: 

a = 1- aKc ■ 
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(12) 



Equation 4.21 in (von Caemmerer, 2000) 



CO, concentration in the bundle sheath Cm 



;+K n 1 + 



Obs)A + R u . 
K 0 



Vc, 



A + Rg 



(13) 



Equation 4.1 1 in (von Caemmerer, 2000) 



Q 

9bs 

K c 

K 0 

K P 

0 BS 



C0 2 partial pressure in M (at the site of PEP carboxylation) C M = Q 
C0 2 concentration in the intercellular spaces as calculated by the IRGA 
Bundle sheath conductance to C0 2 
RuBisCO Michaelis-Menten constant for C0 2 
RuBisCO Michaelis-Menten constant for 0 2 
PEPC Michaelis-Menten constant 

0 2 mol fraction in the bundle sheath cells (in air at equilibrium) o BS = 0^ + 
0 2 partial pressure in the mesophyll cells (in air at equilibrium) 



aA 



0.047g B; 



(14) 



(ibar 
Hbar 

0.005 mol m 2 s~ 1 

650 (ibar (von Caemmerer, 2000) 
450000 (ibar (von Caemmerer, 2000) 
80 |ibar (von Caemmerer, 2000) 

(imol mor 1 Equation 4.16 in 

(von Caemmerer, 2000) 
210000 (ibar 



Flight 

^Cmax 

Vo/Vc 
V P 

^Pmax 



Y 

r 



Respiration in the light, assumed to equal dark respiration 

Mesophyll non photorespiratory C0 2 production in the light fi M = 0.5 R L ight 



Maximum RuBisCO carboxylation rate 



CruiVpr 



Vq_2T 
C BS 



PEP Carboxylation rate V P 

Cm + K P 

Maximum PEPC carboxylation rate 
Fraction of PSII active in BS cells 

Half of the reciprocal of the RuBisCO specificity 

C0 2 compensation point in absence of dark respiration r* = y'Oes 



(imol nr 2 s~ 1 (von Caemmerer, 2000; 
Kromdijk ef al. , 201 0; Ubierna ef al. , 201 3) 
60 (xmol nr 2 s~ 1 (von Caemmerer, 2000 
(1 5) Equation 4.8 in (von Caemmerer, 2000) 



(1 6) Equation 4.1 7 in (von Caemmerer, 2000) 

120 (imol rrr 2 s~ 1 (von Caemmerer, 2000) 
0.15 (Edwards and Baker, 1993; von 
Caemmerer, 2000; Kromdijk et al. , 2010) 
0.000193 (von Caemmerer, 2000) 

(1 7) Equation 4.9 in (von Caemmerer, 2000) 



dark respiration on each individual plant after the measure- 
ments in the light. 

As shown in Fig. 3, errors in the determination of Fm' sug- 
gest that techniques such as the multiphase flash (Loriaux 
et al, 2013), or initial checks to ensure that the saturating 
pulse is saturating (see Bellasio and Griffiths, 2014b) are 
normally appropriate for this method. However, the use of 
our method is possible without a multiphase flash. Firstly, 
the underestimation of Fm' introduces a systematic error, 
i.e. comparable plants will normally show similar VqIV c 
(see Bellasio et al, 2014), unless the extent of C 4 or C 2 
activity has changed under these conditions. Thus, the 
precision and the resolution of the method, when com- 
paring different phenotypes against a common genetic 
background, are not affected by a consistent underes- 
timation of Fm'. Secondly, to improve accuracy, i.e. 



the capacity of the method to estimate the true V 0 /V c , 
other approaches could: (i) increase the saturating pulse 
intensity; (ii) reduce the distance between light source or 
fibre-optic probe and leaf (in some systems); (iii) decrease 
actinic light intensity (as shown in this study) to maxim- 
ise Y(II); and (iv) C0 2 concentration can be increased, in 
order to maximise YfllJ. 

IRGA recalibration, matching Y(ll), Q, and 
consideration of mesophyll conductance 

As mentioned in the results, a slight effect on stomatal con- 
ductance and Q (under low 0 2 ) could have been caused by 
not recalibrating the IRGA upon switching background 
gas (Bunce, 2002). Although that recalibration could have 
increased Q and g s accuracy (under low 0 2 ), this procedure 
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is liable to introduce operator error and extend the time taken 
for measurements; further, there are theoretical reasons why 
we need not account for these processes while carrying out 
such a simple comparative screen. Firstly, the data used to 
calculate equation 8 are measured by the C0 2 channel of 
the IRGA and the fluorometer, which are both unaffected by 
the background gas (Bunce, 2002). Secondly, the effect of Q 
on A (under low 0 2 ) is, for the greatest part, accounted by 
the feedback on Y(II). Although Q decreases under low 
0 2 , there is a strong feedback between assimilation and 
Y(II), and therefore Y(II) decreases proportionally. In 
fact, the relationship between gross assimilation (or, bet- 
ter, between Y(C0 2 ), which is GA divided by PPFD) and 
Y(II) is strictly linear (Edwards and Baker, 1993; Valentini 
et al, 1995; Martins et al, 2013). In C 4 plants, this lin- 
ear relationship has generally a zero intercept, (Edwards 
and Baker, 1993); therefore, for C 4 plants, there is no need 
for curve fitting and the relationship can be correctly esti- 
mated with a single point. In C 3 systems this relationship 
is still linear but the intercept is, although generally small, 
not zero. The intercept, which is the magnitude of engage - 
mant of alternative sinks, can be estimated by linear curve 
fitting, although several data points are required (Valentini 
et al, 1995; Martins et al, 2013). Using the complete fit- 
ting of the Y( C0 2 )l Y( II) relationship, however, did not 
improve the estimate of V Q /V C (data not shown): the com- 
plete curve fitting correctly estimates the intercept, but the 
datapoints are taken under conditions which differ from 
those under which VJV C is measured. 

Another way to improve the estimate of V Q /V C would 
be to adjust C a under low 0 2 so as to match Y( II) meas- 
ured under ambient 0 2 with Y( II) measured under low 0 2 . 
Alternatively, C a could be manipulated to deliver C; under 
low 0 2 , which matches that under ambient. The advantages 
would be that the measured data would then probably fit 
the predicted C 3 and C 4 models more precisely when Q is 
limiting (see Fig. 2, Tables 3 and 4). However, these oper- 
ations do not improve the capacity to screen between C 3 
and C 4 photosynthesis and the additional manipulations 
increase time and likelihood of errors. We also note that 
such improvements would allow this method to be used to 
calculate the C0 2 concentration at the site of carboxyla- 
tion (C c ) in C 3 plants through equation 1 1 (Table 3), as well 
as mesophyll conductance via equation 10, using C c , and 
the values for assimilation and Q measured under ambient 
conditions. 

Conclusion 

In this paper a simple method, and associated theory, have 
been presented, which allow the determination of both the 
oxygenation (V Q ) and carboxylation (V c ) rate of RuBisCO 
and the rate of photorespiratory C0 2 evolution (F) based 
on gas exchange and variable chlorophyll fluorescence under 
ambient and low 0 2 . This may be of particular interest for 
high throughput screening to identify C 4 mutants lacking a 
fully functional CCM, C 2 variants, or populations of C 3 -C 4 
hybrids (Oakley et al, 2014). 
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